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Abstract. The collision energy dependence of jet tomography is investigated 
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1. Introduction 

We study the energy systematics of jet tomography in nuclear collisions within the 
GLV formalism[n Q] m the y/s — 62 .4 to 200 AGeV range. Jet tomography at 
\fs = 62.4 AGeV data [0] tests the predicted y/s decrease of the QGP density and 
the predicted variation of the gluon/quark jet source. In this paper, we calculate 
the nuclear modification factor, Raa(j>t, V = 0, y/s), for central Au — Au collisions 
at both y/s = 62.4, 200 AGeV for neutral pions. Previous predictions for 62 AGeV 
have been published by Wang[0] and Vitev[EJ. 

We concentrate here on the role of energy loss fluctuations [0 E] on the pre- 
dicted single hadron attenuation pattern in order to gain an estimate of some of the 
systematic theoretical errors in the jet tomographic technique. To isolate the role 
of fluctuations we neglect kx smearing, the Cronin enhancement, gluon and quark 
shadowing, and nonperturbative baryon dynamical contributions that strongly dis- 
tort the hadron spectra below pr < 4 — 5 GeV[Q. So our results are valid for 7r° 
spectra for pr > 4 — 5 GeV. 

We test the influence of the shape of the energy loss fraction spectrum, P(e, e) 
about the mean energy loss fraction, e, which was first pointed out as important in 
[ |Hj . This sensitivity to the shape increases as the high pt slopes increase at lower 
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v/i. 

2. Calculation of Spectra and R AA 

The neutral pion cross section in pp can be calculated using collinear factorized 
pQCD. In a dense QCD medium the induced radiative energy loss reduces the 
initial px of the jet parton by a fraction e before hadronization. In this paper we 
calculate the tt° inclusive spectrum as the following. 

E h—# =A Z^ dx 1 dx 2 def a/A (x 1 ,Q )f b/A (x2,Q ) -j — P{e, e) 

abed 

(1) 

where z* = j^The inclusive number distribution is calculated by mutiplying 
this invariant distribution with the the Glauber geometric binary collisions factor, 
TaaQ))- The functions f a /A ancl Dh/c are the conventional MRS D- distribution 
function and KKP fragmentation function, respectively. Gluon number fluctuations 
are taken into account using the distribution P(e, e), where e is interpreted as the 
average fractional energy loss and is proportional the local gluon rapidity density. 
We explore two simplified forms of fluctuation distributions to assess some of the 
systematic uncertainties in the predicted quenching. One is a "uniform" model 
that essentially reproduces the truncated Poisson of [EJ 0. The second is called 
"squeezed" because it accumulates strength near the e « 1 opaque limit. This 
distribution is considered to take into account the alternative branching form of 
implementing gluon fluctuations. 

The "uniform" distribution takes the form, 



fl(0<e<2£) 



if < e < 0.5 



while the "squeezed" distributon is the following. 

f «(o«l<*-) if o < e < 0.5 



2(1- 



The factors of K and Q 2 are fit to pp data at the requisite COM energies, e evolves 
with y/s according to the the multiplicity evolution. We assume that e c (y/s) — 



Tf~ ( Tn° (wo) /dy ) e a (200) where c is the parton type and C c / g are the QCD Casimirs. 
Thus, the free parameter is e g (200) which is set to fit PHENIX [0 tt data at 200 
AGeV. 

Once the spectra for the A — A and p — p reactions at the requisite energies 
have been calculated, the nuclear modification factor Raa is just the ratio between 
them. This gives a range of values for each distribution ( "uniform" and "squeezed" ) 
determined by the errors of the data. For the "squeezed" distribution, the deter- 
mined range is 0.65 < e s (200) < 0.76 while the "uniform" distribution can be fit 
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to Raa data with the range 0.70 < e g (200) < 0.80. The "squeezed" distribution 
needs a lower average opacity as the distribution itself is biased towards e — > 1 when 
e > 0.5. 
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Fig. 1. Plots of Mepredicted bands of Raa{pt) at ^fl = 62.4'^GeV for both tie 
"squeezed" and "uniform" distributions. The figure on the left shows the predicted 
Raa for the "squeezed" distribution while the figure on the right shows Raa for 
the "uniform" distribution. 



The predicted Raa at 62.4 AGeV is obtained using the multiplicity systematics 
from PHOBOS [HHj that suggests the value dN g /dy{^ = 62.4) « 650 - 770. We 
expect from hadron gluon duality arguments that this multiplicity is reduced from 
about 1000 at 200 AGeV. Therefore, we expect that e(62.4) « (0.65 - 0.77)e(200). 
The bands found by fitting to 200 AGeV PHENIX data can now be extrapolated to 
further bands at 62.4 AGeV. The predicted Raa(pt) for both distributions can be 
seen in Fig. ^ Using our multiplicity extrapolations we find that for the "squeezed" 
distribution and dN g /dy(62A) w 650(770), 0.42(0.50) < e g (62.4) < 0.49(0.59). 
Similarly, for the "uniform" distribution and dN g /dy(62A) « 650(770), 0.46(0.54) < 
e ff (62.4) < 0.52(0.62). One of the things to note is that the "uniform" distribution 
prediction is significantly flatter over px than its "squeezed" counterpart. This is 
because once e > \ the two distributions treat quenching very differently. The 
"uniform" distribution saturates to a uniform distribution over < e < 1 while the 
"squeezed" distribution piles up closer to e — > 1 and causing a larger variation in 
the quench. The calculations for the nuclear modification at s/s — 62.4 AGeV are 
consistent with Vitev Any deviations between the spectra can be attributed to 
the inclusion of models for Cronin interactions as well as initial parton kx smearing 
in [0 which are not included in the current calculations. 

Note that the predicted Raa(pt, V~s = 62.4) (Fig. ^) have a negative pr slope 
compared to the generally flat Raa at 200 AGeV. This higher slope is caused by 
the earlier set in of the kinematic limits of the problem at lower energies. The 
"kinematic suppression" can be more robustly seen by calculating the observable 
R s (s) — j^^200) seen * n ^ ^ s curves have a distinct downwards slope due 
to the increasing power of the initial parton distributions. R s is perhaps a better 
observable to use than Raa to observe the energy dependence of jet quenching as 
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Solid Lines indicate dN/dy = 771) 




^-s Dashed Lines indicate dN /dv = 650 






1 3 denote Squeezed Distribution 




2,4 denote Uniform Distribution 
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Fig. 2. The ratio for R AA at 62.4 AGeV to that at 200 AGeV. Curves for "squeezed" 
distributions are generally higher than curves for "uniform" distributions. 

the uncertainty in the multiplicity extrapolations get canceled in the ratio. R s also 
able to differentiate between the two types of fluctuation distributions (see Fig. |5J). 
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